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ABSTRACT 

General characteristics were determined for a gas- coded reactor designed to 
operate at 1 MW for 25 000 hours, with a coolant outlet temperature of 1650 S, Design 
curves relating core size to fuel concent ration are presented as a function of radial 
power flattening. 
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A 1-MEGAWATT REACTOR DESIGN FOR BRAYTON-CYCLE 
SPACE POWER APPLICATION 
II -NEUTRONICS DESIGN 

by Charles L. Whitmarsh, Jr. , and PaulT. Kerwin 
Lewis Research Center 

SUMMARY 

Neutronics design calculations were performed for a 1 -megawatt thermal gas -cooled 
reactor required to operate for 25 000 hours with a coolant outlet temperature of 1650 K. 
The reactor model is a cylindrical core (length = diameter) contained in an 0.95- 
centimeter -thick tungsten pressure vessel. Core structural materials are tantalum, 
tungsten, and their alloys. The fuel is either enriched uranium dioxide or uranium 
nitride. Control is achieved by axial movement of beryllium oxide radial reflectors. 

The control range is a nominal 10 percent reactivity with a maximum multiplication 
constant of 1. 05. 

Design curves relating minimum core size to uranium concentration are generated 
for reactors reflected axially by a 1.27-centimeter-thick tungsten end plate and radially 
by a 7.62-centimeter-thick beryllium oxide reflector. These design curves include 
radial peak- to average-power ratios of 1.05, 1. 10, and 1. 35. In addition parametric 
data are presented such that design curves for other conditions can be calculated. 

The radial power flattening (achieved by variation of fuel enrichment in annular 
zones) sufficiently enhanced reflector worth so that control is feasible inside an adjacent 
4 it shield. Spatial and energy flux distributions calculated for a representative configu- 
ration indicate that the median fission energy and prompt- neutron lifetime in the core 
are 0. 35 MeV and 0. 3 microsecond, respectively. 


INTRODUCTION 

Nuclear- reactor -powered Bray ton- cycle systems are of interest for space applica- 
tions requiring electric power in the range of tens of kilowatts to megawatts (refs. 1 and 
2). The Brayton system has the advantages of a noncorrosive working fluid, high cycle 



efficiency , and a relatively well- developed component technology; but it also has such 
disadvantages as large piping and components to be shielded and large radiator area. 

The application of this reactor system at high powers is generally restricted by radiator 
area. 

Although either a liquid-metal- cooled or a gas -cooled reactor is useful as a heat 
source for Bray ton systems, only gas -cooled reactors are considered herein. Such 
design goals as high- temperature operation and small core size have led to typical 
reactor characteristics of highly enriched, highly dense fuel, refractory-metal struc- 
tural components, and reflector control systems (ref, 3), As a result the cores have 
high-energy flux spectrums and are categorized as "fast" reactors. 

This study is of gas- cooled reactors for use in an M advanced-state-of-the-art M 
Brayt on- cycle space power system. The objective is a preliminary reactor design that 
would produce 1 megawatt of thermal power for 25 000 hours with a turbine inlet temper- 
ature of 1650 K and that would retain all fission gases in the core, A single -loop system 
was selected so that the reactor coolant outlet temperature could be assumed to be equal 
to the turbine inlet temperature. The reactor core is composed of either fuel-pin bun- 
dles or a solid matrix fuel penetrated by coolant passages. Since reactor size with its 
associated shield weight penalty is of major importance to space power systems, design 
effort was concentrated on minimization of core size. 

The design study is presented in two reports: I - THERMAL ANALYSIS AND CORE 
DESIGN, and n - NEUTRONICS DESIGN, Part I contains the heat-transfer and fluid-flow 
calculations required to define the internal core configuration and to show the effect of 
various design parameters on core size. Part n defines the general reactor configura- 
tion, the control system, and the materials of construction and presents core sizing data. 

More specifically, this report, part H, NEUTRONICS DESIGN, shows the relation 
between core size and required fuel content for a range of nuclear design variables. The 
primary variables are reflector worth and radial power flattening. In addition, the 
neutron-flux spatial and energy distributions, the neutron lifetime, the control system, 
and the shield effects which are typical of this general class of reactors are presented. 


REACTOR CONCEPT 

A cycle analysis (ref, 4} for a one-loop Brayton-cycle system, which uses a 1- 
megawatt nuclear reactor heat source to supply 1650 K gas to a turbine, resulted in the 
recommended conditions shown in figure 1. For the cycle analysis the reactor was 
treated as merely a heat source with only the inlet and outlet conditions specified- For 
the neutronics analysis a cylindrical core with a length to diameter ratio (L/D) of 1 was 
selected because it provides (1) small critical core size, (2) good heat-transfer charac- 
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Figure 1. - Schematic drawing of Brayton-cyde system, Workifig fluid, 40 percent He - 60 
percent Xe gas mixture (molecular weight - SO kg/kg-mde); mass flow rate, S.5 kilograms 
per second. 


teristics, (3) ease of fabrication, and (4) application to either shadow or 4ir shielding 
systems. The core could consist of either (1) bundles of fueled tubes with the coolant 
gas external to and flowing parallel to the tubes (pin core) or (2) a number of hexagonal 
fuel elements fitted together and penetrated by holes through which the coolant gas flows 
(matrix core). These concepts are illustrated in figures 2(a) and (b), respectively. 

Two fuel arrangements in the core were considered, namely, uniform distribution 
of fully enriched fuel (93. 2 percent U c ) and nonuniform distribution in which the core 
was divided into six radial zones with fuel enrichment varied among zones to achieve a 
flat power distribution (ref. 5). Space for the collection of fission gases is provided in 
the fuel pins but was not needed in matrix fuel elements. The core is contained in a 
pressure vessel and has end plates at both ends. Control is effected by a split radial 
reflector which forms a sheath around the core. Leakage of neutrons (and therefore 
reactivity control) can be regulated by axial movement of all or part of the reflector 
(fig. 3). 

For a small core, an external control system is preferred, for example, reflector 
control. However, the particular method of operation (leakage control by axial move- 
ment) was chosen somewhat arbitrarily. A literature search (refs. 6 to 8) indicated that 

(1) radial reflectors could give sufficient reactivity changes by regulating core leakage, 

(2) end plates were effective as axial reflectors, (3) additional axial reflectors saved 
little in overall length (core plus reflector), and (4) heavy reflector materials were 
almost as effective per unit of thickness as light (moderating) materials. 

The choice of materials is primarily dependent on the high-temperature (>1900 K) 
environment in the reactor. Selection of structural materials was thus limited to 
tungsten (W), tungsten - 25 percent rhenium alloy (W-25Re), and a tantalum-base alloy 
(ASTAR-811C) (ref. 9). Uranium dioxide (UOg) (ref. 10), uranium nitride (UN) (ref. 11), 
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(a) Fully reflected end-of- (b) Operating control (d Alt reflectors open, 

life configuration. reflectors open, beginning- shutdown configuration. 

of-life fhot, clean, critical) 
con fig oration. 

Figure 3, - Radial reflector configuration. 



Figure 4 . - Schematic drawing of gas -cooled, 
reflector-controlled, fast- sped rum reactor. 


and a matrix of 60 volume percent UO 2 in W {ref. 12) were considered for reactor fuels. 
These materials have good high- temperature properties and high uranium concentration. 
In addition to having a high-temperature capability, the reflector material must also 
have a small neutron-absorption cross section for fast neutrons. Beryllium oxide (BeO) 
was selected over several applicable materials, such as molybdenum (Mo), niobium (Nb), 
or alumina (AlgOg) because of lower weight for equivalent reactivity worth. The reactor 
concept is shown in figure 4. 


ANALYTICAL PROCEDURE 
Calculationa! Model 

Based on the concept developed in the preceding section, the geometric model shown 
in figure 5 was used in the reactor calculations. Only a quarter of a vertical cross- 
section is shown because the reactor configuration is symmetrical about the axial 
centerline and reactor midplane. The same model was applicable to either the pin or 


Pressure vessel 


Coolant plenum 


End plate 


Fission -gas chamber 


-Axial centerline 


Core 


Radial midplane - 


Reflector 


“Gap (0. 635 cm) 


Figure 5. - Neutronfcs calculation model. 
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the matrix core calculations because the core region was assumed to be homogenized. 
This assumption permits use of one -dimensional calculation methods and is an accurate 
assumption when neutron mean free paths are large compared with region dimensions. 
The complete absence of any in-core moderating materials indicates that all designs will 
have a high-energy neutron spectrum, which in turn means that neutron mean free paths 
will be relatively large (>2 cm). 


Criticality Calculations 

The anticipated small core size (diameter of 30 to 40 cm) and fast-neutron spectrum 
led to the use of the multigroup transport codes DTF (ref. 13) and TDSN (ref. 14). A 
set of 16-group spectra-averaged cross sections was generated with the GAM-n (ref. 15) 
and TEMPEST- n (ref. 16) codes. Cross sections and code options were verified by cal- 
culations of critical experiments with similar materials and configurations (see 
appendix A) . 


Reactivity Calculations 

Design changes, for example, a change in reflector thickness, can be represented 
by changes in reactivity. To be useful as parametric design data, these reactivity effects 
should be both independent of each other when more than one design change is made on 
the same system and applicable to a noncritical system. The reactivity equations 
utilized in this report are as follows: 




k 2 ~ k l 
k 2 k l 


x 100 


N 


A Ptotal = y>i 


i=l 


( 1 ) 


( 2 ) 


where Ap^ is the i th component of the total percent reactivity change (or reac- 
tivity worth), kj is the multiplication constant of the initial reactor, and kg is the 
multiplication constant of the perturbed reactor as derived from perturbation theory. 
Thus, the results are subject to the same assumptions and limitations as perturbation 
theory (ref. 17). The more important limitation is that the flux energy spectrum is 
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assumed to be constant from the initial to the perturbed state. Accuracy, therefore, 
diminishes with increasing deviation from criticality (k = 1). 


RESULTS 

Criticality 

Because most nonfissile reactor materials have similar fast-neutron cross sections, 
criticality calculations are somewhat insensitive to small changes in composition. Thus, 
parametric calculations based on cores with nominal compositions (void volume fractions 
f of 0. 15 and cladding- to fuel- volume-fraction ratios f c /ff of 0. 1 to 1. 0) may be used 
for general core sizing. The relation between U 235 density in the core and the multipli- 
cation constant for cylindrical, "bare” cores with L/D = 1 and diameters of 29 to 42.7 
centimeters is shown in figure 6 . Bare core is used here to denote an unreflected core 
surrounded radially by a 0. 95-centimeter-thick tungsten pressure vessel. All calcula- 
tions were based on fully enriched (93.2 percent U 23 ^) UO 2 or UN fuel with either 
ASTAR-811C or W-25Re as the fuel cladding material. 

Although pressure vessel thickness generally is dependent on diameter, use of a 
constant thickness for neutronics calculations appears justified, hi a fast-spectrum 
reactor, the pressure vessel acts more as a reflector than an absorber, and since its 



Figure 6. - Multiplication constant as function of 
uranium-SS concentration for tin reflected, cy~ 
liftdrlcai, fast-spectrum reactor* Pressure vessel, 

0, 95-cefttimeter tungsten- length to diameter ratio, L 
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thickness is small compared with the actual reflector, the multiplication constant is 
insensitive to small changes in pressure vessel thickness. A value of about 1 percent 
Ak gain per centimeter increase in thickness was calculated for a tungsten pressure 
vessel. 


Control System 

Reactivity requirements. - The reactivity required for reactor operation can be 
classified into four categories: (1) fuelburnup, (2) Doppler effect, (3) temperature 
effect, and (4) shutdown margin. 

Fuel burnup was calculated from the design specifications of 25 000 hours operation 
at 1 megawatt as follows: 


Burnup = Power x Time 
or 

B™ = !0 6 watts x 3. lxlO 10 fissions x 22 P g 1)235 x — 1 

watt- sec mole 6. 023X10 23 atoms U 235 

x x 9X10* 7 sec (3) 

fission 


B t = 1100 g u 235 

235 

For the critical masses anticipated in this study, 1. 1 kg of U represents about 0. 5 to 
1 percent of the inventory of large and small cores, respectively. This fuel loss repre- 
sents a reactivity change of about 0. 5 percent (see fig. 6). Because of the smallness of 
the reactivity change, a more detailed calculation including fission-product buildup was 
considered unnecessary. Thus, a nominal 1- percent reactivity allowance was provided 
for fuel burnup. 

Calculations of Doppler coefficient were not attempted because accurate values are 
difficult to obtain for fast, highly enriched reactors and are therefore considered beyond 
the scope of a parametric study such as this. The Doppler coefficient probably is nega- 
tive but sufficiently small (ref. 18) that zero was assumed. Approximate calculations 
indicated values near -0. 1 percent reactivity for small cores and -0.5 percent reactivity 
for large cores. 

The temperature effect represents the change of reactivity from the cold-clean con- 
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diiion to full-power operation resulting from material density and core geometry changes. 
Criticality calculations on a representative core (D = 32. 6 cm, f^ = 0, 72, = 0.073) 

indicated a loss of ~1, 8 percent reactivity because of these changes. These calculations 
were based on estimated average operating temperatures of 1720 K for the fuel, 1520 K 
for the cladding and pressure vessel, and 1250 K for the reflector. Homogenized cores 
were used in the calculation and therefore detailed geometric effects within the core 
were not considered. However, expansion of the reflector and change in core length were 
included in the calculation. 

The amount of suberitical reactivity needed for shutdown margin is primarily a 
matter of the safety and hazards analysis and is dependent on the application. One 
criterion which seems reasonable for this study is that with half the control system in 
the "in ,, position, the cold -clean core should be suberitical. Thus, if a 2. 2 percent 
reactivity contingency is added to the temperature and burnup reactivity requirements, a 
shutdown margin of about 5 percent reactivity is required. It follows that the design 
multiplication factor will be 1.05 and the total control range of the reflectors must be 
greater than 10 percent reactivity (table I). 


TABLE I. - REACTIVITY CONTROL 
REQUIREMENTS 


Category 

Nomina! 

reactivity 

needed, 

percent 

Fuel burnup 

1.0 

Doppler 


Temperature delect 

1.8 

Contingency 

2.2 

Shutdown 

5,0 

Total 

10.0 


Reflector worth . - For a given core diameter, reactivity worth of the reflector in- 
creases with increasing reflector thickness, as shown in figure 7. However, the rate of 
increase lessens as reflector thickness increases. For a fixed reflector thickness, 
reactivity worth decreases with increasing core diameter. These data represent 
criticality calculations for uniformly fueled cores surrounded by a 0. 9 5- centimeter - 
thick tungsten pressure vessel both with and without a radial reflector. Also shown in 
figure 7 is a relative increase in reflector worth of 0. 5 (about 6 percent reactivity) for a 
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Figure 7. - Reactivity worth of radial BeO re- 
flectors surrounding cylindrical, fast- 
sped rum reader core. Fuel, OOg; cladding, 
ASTAR-mC; pressure vessel, 0.95- 
centi meter tungsten; length to diameter 
ratio, 1. 


32. 6-centimeter-diameter core when fuel enrichment is varied in an nul a r zones to obtain 
a radial peak- to- average power ratio (P p /P) R of 1. 05 (case H, ref. 5). 

Thus, the minimal control requirement of about 10 percent reactivity should be 
obtainable for uniformly fueled cores less than 34 centimeters in diameter with a radial 
BeO reflector greater than 5 centimeters thick. Larger cores would require thicker 
reflectors. Operating control of 1. 0 percent reactivity required to compensate for fuel 
burnup over the core life could be regulated by moving the outer 2 centimeters of the 
reflectors of an unzoned core and the outer centimeter for a zoned core. 

The reactivity worth of a radial reflector as a function of eore diameter is shown in 
figure 8. These results were calculated for uniformly fueled, fully enriched cores with 
different compositions, although all had a 7. 62- centimeter BeO radial reflector. 

The effect of axial reflection was simulated by varying the thickness of the tungsten 
end plates. Symmetric configurations with the same thickness of tungsten at each end of 
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figure 8 . * Reactivity worth of 7, 62 -centimeter BeO radial 
reflector surrounding cylindrical, uniformly fueled, fast 
reactor core with length to diameter ratio of 1 . 



End-plate thickness, cm 


Figure 9. - Reactivity worth of tungsten end 
plates in BeO reflected, fast-spectrum, 
cylindrical reactor core. End plate, 80 per- 
cent W - 2 fl percent void; fuel, UO 2 ; clad- 
ding, A5TAR-8I1C; pressure vessel, 0.95- 
centimeter tungsten; length to diameter 
ratio, 4 7. 62 -centimeter 6eO radial re- 
flector. 


the core were considered, and the tungsten density was reduced 20 percent to account for 
coolant passages. The shape of the reactivity-worth curve was similar to the one for 
the BeO radial reflectors although the worth per unit thickness was less (fig. 9). The 
qualitative results of radial and axial reflectors should be the same, that is, decreased 
worth with increased core size at constant length to diameter ratio. However, calcula- 
tions indicated axial reflector worths to be relatively insensitive to such changes. For 
a core-diamter change from 30 to 42 centimeters, the change in worth of the 1. 27- 
centimeter tungsten axial reflector was less than 0.6 percent reactivity. The effect on 
axial reflector worth of c han ging radial reflector thickness from 7.62 centimeters to 
5. 08 centimeters on a 32. 6 -centimeter -diameter core was calculated to be of similar 
magnitude. Therefore, a single curve (fig. 9) was used to determine axial reflector 
worth for all core configurations considered in this study. 

The final determination of end- plate thickness will probably depend on fabrication 
considerations rather than reactivity. With this in mind, a value of 1.27 centimeters 
was selected in this parametric study for the reference design calculations. 

A comparison of radial reflector worths of different materials showed that, for a 
reflector thickness of 6.35 centimeters, the moderating material BeO is worth 3 to 4 
percent more reactivity than the nonmoderating materials Mo or Nb (table n). Part of 
this increase was due to the greater n-2n reaction in Be, which indicates that the reflec- 
tion properties per unit thickness of these materials are nearly equal. 

Shield effect. - Leakage- controlled fast reactors are inherently sensitive to any 
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TABLE n. - EFFECT OF REFLECTOR TABLE HI. - EFFECT OF INFINITELY THICK TUNGSTEN 

MATERIAL ON REACTIVITY WORTH SHELL (4 ir SHIELD) ON UNZONED FAST -SPECTRUM 

OF 6. 35-CENTIMETER RADIAL REACTOR REFLECTED WITH 7.62 


REFLECTOR ON SMALL, 


CENTIMETERS OF BeO 


ZONED -CORE, FAST- 
SPECTRUM 


REACTOR 


Reflector 

Reactivity worth 

material 

Percent ; 

Relative 

BeO 

19 

1,0 

Mo 

16 

.84 

Nb 

15 

,79 


Distance from oater 
1 surface of 

reflector, 
cm 

Control range with 
7,62«cm BeO shell, 
percent reactivity 

Center- to average- 
power ratio 

(No shield) 

22 

1,41 

30,5 

19 

1.38 

15,2 

17 

1.37 

0 

7.6 

(a) 


a Not calculated but should be about 1.37. 


material external to the reactor system since this material acts as a neutron reflector. 
The effect of shield material on the reflector control range was therefore investigated 
because application of these reactor systems would require shielding. Since specific 
shielding requirements are unknown, this investigation was both general and brief. 

Results presented in table HI indicate the reactivity control range of a 7. 62- 
centimeter BeO reflector when an infinitely thick tungsten 47? shield is placed at various 
distances from the outer edge of the reflector. The two extreme positions, no shield 
and an adjacent shield, show a difference of 14 percent reactivity in the control range of 
a 7. 62-centimeter BeO reflector. Other effects noted were an increase of multiplication 
constant k and a reduction of the center- to average-power ratio from 1. 41 to 1. 37. 

These results are based on a spherical reactor with an 18. 5- centimeter- diameter 
core surrounded by a 0. 95-centimeter-thick tungsten pressure vessel and a 7. 62- 
centimeter BeO reflector. Therefore, these results are only qualitative, and the trends 
are applicable to cylindrical systems with radial reflectors. 

Since the control range of an unzoned spherical core with an adjacent shield was less 
than the nominal 10 percent reactivity required in this study, a more specific calculation 
was performed. The reflector control range of a reference-design core (33, 5-cm-diam 
cylindrical core with L/D = 1, radially zoned to (Pp/p) R = 1.05, a 7.62-cm-thick BeO 
radial reflector, and k = 1.05) was calculated to be 8. 4 percent reactivity with an 
infinite tungsten shield adjacent to the BeO radial reflector. This value is sufficiently 
close to the 10 percent reactivity nominal control requirement to indicate design feasi- 
bility. However, other factors are involved which would require a redesign of the core 
for operation inside a shield. The reference core power ratio (P n /P) R increased to 
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1. 12 in the peripheral fuel zone and the design k increased to 1. 089, These two effects 
on core size are opposing, and the net result is unknown without additional calculations. 
For this study the additional design effort was considered unwarranted, and it should 
suffice to say that operation with an adjacent shield appears feasible. 


Power Distribution 

The radial power distribution calculated for a uniformly fueled, fully enriched reac- 
tor with a BeO radial reflector is presented in figure 10. The local- to average-power 
ratio, which peaks at the core center, decreased from 1.42 to 1.28 as the radial reflec- 
tor thickness was increased from 5. 08 to 7. 62 centimeters. These reflector thicknesses 
represent probable initial and end-of-life positions. There was no indication of power 
peaking at the core -reflector interface. This design problem is often encountered when 
a fast core is adjacent to a moderating medium since high-energy neutrons that escape 
from the core and return at lower energies are immediately absorbed by the dense fuel 
material. Because of anticipated materials problems related to fission-gas retention in 
core regions with relatively high power output (ref. 4), power flattening by variation of 
fuel enrichment in radial core zones was considered. The neutronics aspects of power 
flattening were reported previously (ref. 5) and are summarized as follows. 

Power flattening by variation of fuel enrichment perturbs the reactivity of a fully 
enriched, uniformly fueled reactor in two ways: 

(1) Reactivity is reduced because fuel must be removed. 



Figure 10. - Radial power proflies in zoned and unzoned fast -spectrum reactors with 8e0 radial re- 
flectors. 
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TABLE IV. - REACTIVITY COST TO RADIAL POWDER FLATTEN 


CYLINDRICAL, REFLECTOR-CONTROLLED, FAST -SPECTRUM 
REACTOR WITH 32. 6 -CENTIMETER DIAMETER AND LENGTH 


[Radial BeO reflector thickness, 5. 08 cm. ] 


Radial peak- to 

Reactivity 

average-power 



■ 


r^tin 

Loss from XT 

Gain from increased 

Net 

/ i~\ 

removal j 

reflector worth 

loss. 

\ r ,*L 

percent 



percent 



Percent 

Relative 


1.05 

12 

6 

1,5 

6 

1.10 

6 

2 

L15 

3 


(2) Reactivity is Increased because the worth of the radial reflector is increased. 
Reactivity values for these items based on a representative core design are listed in 
table IV. This design assumes that shutdown control will be achieved by a 5. 08- 
centimeter BeO radial reflector and operating control by an additional 2. 54 centimeters 
of BeO. The peak- to average-power ratios of 1. 05 and 1. 10 were estimated to be lower 
limits that could be reasonably obtained by fuel zoning. Thus, the reactivity cost to zone 
to a (Pp/P)^ of 1* 05 was 6 percent and to a (P^/P)^ of 1- 1 was 3 percent. 

These data are indicative of the initial (beginning of life) core configuration only, 
with the shutdown reflector closed and the operating control reflector open. The result- 
ing distortion of the radial power shape caused by reflector movement to compensate for 
fuel burnup over the core life was calculated for an operating reflector thickness of 2. 54 
centimeters. This nominal thickness was selected in the initial stage of the study and 
therefore does not coincide with the outer 1 centimeter indicated in the section R eflec tor 
worth. However, any effects of this difference in reflector thickness on the final core 
design are considered to be small. 

The data in table V (and fig. 10) indicate that the (P p /P) R in the two peripheral zones 
exceeds 1. 05 in the fully reflected (7. 62 cm) configuration. However, for some design 
calculations, such as fission -gas production, the quantity of interest is the time -averaged 
(or "effective”) value of (P p/P) R - H the operating reflector is assumed to be withdrawn 
at the beginning of core life and inserted at the end of core life and its movement is 
assumed to be linear with fuel burnup, the effective value for (P^/P) R is then about t. 05 
and occurs in the outer two zones. 

An alternate power -flattening method is to use a fuel distribution which would pro- 
vide a (P p/P) R of 1. 05 for the end-of-life configuration (7. 62-cm BeO radial reflector). 
Data presented in reference 5 show that, although the reactivity penalty would be less 
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TABLE V. - DISTORTION OF RADIAL POWER SHAPE IN ZONED -CORE, 
FAST -SPECTRUM REACTOR CAUSED BY CHANGING 


REFLECTOR THICKNESS 


Radial &one 

Enrichment, 

percentage of 
jj235 

Initial reflector 
configuration, 
$.08 centimeters 
BeO 

Bnd* 0 f~iile 
reflector 
configuration 3 ', 
7.62 centimeters 
BeO 

Radial peak 

- to average- 

-power ratio, (p /PI , 
v P 'r 

i (central) 

SI. 9 

1 . 

05 

1,00 

2 

80,2 



1.00 

3 

65.8 



1.02 

4 

71.8 



1,04 

5 

81,0 



1.06 

6 (peripheral) 

93.2 

1 . 

03 

1.08 


a Zero fuel depletion was considered in the calculation. An estimated value 
of 1.09 would occur in the peripheral zone if 25 000-MW-hr buraup were 
considered. 


(3 percent compared with 6 percent), power distortion over the core lifetime would be 
greater. For example, the central zone would have an initial (P /P) R of 1. 10 and a 
final value of 1.05, with the effective (or average) ratio being greater than 1. 05. There- 
fore, this method was judged to be less desirable. 

The axial power distribution was calculated to be of a shape that could be approxi- 
mated by cos 7 rZ/Le where Le/L, the ratio of extrapolated length to actual length, was 
dependent on the end -plate thickness {table VI). For a thickness range of zero to 7. 62 
centimeters, Le/L varied from 1. 15 to 1. 56 and the axial peak- to average-power ratio 

(P /p)™ varied from 1. 37 to 1. 19. For the selected thickness of 1. 27 centimeters, the 
P & ^ 

corresponding values were 1. 31 for (P /P) z am * 1- 27 for Le/L. Although the symmet- 
ric cosine function was used throughout this study to represent the axial power distri- 
bution, it was recognized that from the standpoint of heat transfer this was not the best 
shape. A nonsym metrical distribution with its maximum near the core inlet would in- 
crease the average temperature difference between fuel and coolant. This heat-transfer 
benefit associated with the axial fuel zoning was considered to be offset by the calcula- 
tional and fabrication complexities. Therefore, no axial power flattening was attempted. 
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TABLE VI. - AXIAL POWER DISTRIBUTION* IN 
SMALL, FAST -SPECTRUM REACTOR 


Case 

Tungsten 
end -plate 
thickness 11 , 
cm 

Axial peak- to 
average-power 
ratio, 

(' rfih 

Ratio of 
extrapolated 
length to active 
fuel length, 

Le/L c 

1 

0 

1.37 

1,15 

d 2 

1.27 

1,31 

1,27 

3 

2.54 

1.27 

1.35 

4 

5,08 

1.22 

1,48 

5 

7.62 

1. 19 

1.56 


^ased on one -dimensional calculations with radial 
reflection approximated by reflector savings* 
b Axial symmetry with end plate located at each end 
of core. 

Calculated value which fits the power distribution 
to a chopped cosine shape, uZ/L e. 

^Checked with two-dimensional calculation with 
7 ,62 -cm BeO radial reflector. 


Design Correlation 

If the assumption is made that the various reactivity effects considered herein are 
both independent of each other and of the multiplication level of the core, these effects 
can be combined in a general design correlation of bare-core size to the design param- 
eters (e. g. , end -plate thickness, radial reflector thickness, and radial fuel zoning to 
attain power flattening). Thus, 

Ap to tai * 0 + 6) V + y + M (5) 


where 

j8 reactivity worth of 7. 62-cm BeO radial reflector on uniformly fueled, fully enriched 
core (fig. 8) 

6 change in reactivity worth for BeO radial reflector other than 7. 62 cm (fig. 7) 

tj relative increase in reflector worth caused by radial fuel zoning (table IV) 

y reactivity worth of W end plates (fig. 9) 

H reactivity effect due to fuel removal caused by radial fuel zoning (table IV) 
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Then, based on a required multiplication factor of X. 05, 


Ap total “ 


1.05 - kj^ 
1.05 


( 6 ) 


where is the bare-core multiplication factor (fig. 6). Also included in this formula- 
tion is the assumption that certain power -flattening effects are independent of core diam- 
eter (e. g. , the relative increase of reflector worth and the reactivity lost by fuel 
removal). 

Equation (6), then, is a correlation relating bare-core size to required fuel concen- 
tration for a particular set of design conditions. These data are presented in figure 11 
for the recommended conditions of a 7. 62-centimeter BeO radial reflector and 1. 27- 
centimeter W end plates and for three degrees of radial power flattening ((P /P)-^ = 1. 35, 
1. 10, and 1. 05). Figure 11 represents the relation of fuel concentration to minimum 



Figure II. - Criticality design curves for fast- 
spectrum, gas -cooled, ref lector -control led 
reactor with radial peak- to average -power 
ratios of 1,G5, 1. 10 r and 1. 35. Cylindrical 
core; fuel, enriched uranium- cladding, 
W-25Re; pressure vessel, 0.95-centimeter 
tungsten; length to diameter ratio, 1 ; axial 
reflector, 1. 27-centimeter tungsten; radial 
reflector, 7. 62-centimeter BeO; multiplica- 
tion constant, k, 1 05, 
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core size for a gas-cooled reactor designed to operate at 1 megawatt for 25 000 hours 
with a coolant outlet temperature of 1650 K. Such a reactor would have the same reflec- 
tor and pressure vessel configuration described in the preceding paragraph. It should be 
noted that the data in figure 11 represent the optimization (minimum core size) of a 
reactor (core plus control system) to perform a specific job. Thus, the correlation 
includes the interrelating effects of the reflector, core, and fuel. 

Jff curves of uranium concentration against core diameter for other combinations of 
reflector thicknesses should be desired, these could be generated by substituting the 
appropriate reactivities from equation (5) into equation (6) to obtain the required bare- 
core multiplication factor from which the corresponding fuel density is obtained 
from figure 6. However, application of the parametric data to reactors with other than 
7. 6 2 -centimeter BeO reflectors (6 # 0) is probably limited to the unzoned core configura- 
tion (no radial power flattening). This procedure also requires the assumption that the 
increase in reflector worth per unit thickness is independent of core diameter. This is 
not rigorous, but the slope of the reactivity- worth curve between reflector thicknesses 
of 5 and 7.6 centimeters appears to be relatively insensitive to diameter change (fig, 7). 

When using equations (5) and (6) to construct design curves, the term (jS + 6)?? must 
be greater than 10 percent to ensure shutdown margin, and the fuel damage from burnup 
must be evaluated. For the core -size range in figure 11, maximum burnup values for 
UOg fuel were calculated to be 4. 7x10^, 3. 6X10^®, and 3. 2xio^ fissions per cubic 
centimeter for the radial power ratios of 1. 35, 1.1, and 1. 05, respectively; and for UN 
fuel to be 6. 1x10^, 4. 6x10^, and 4. 2x10^ fissions per cubic centimeter. 


Core Characteristics 

Flux spectra . - Flux spectra, calculated at various radial positions in a representa- 
tive zoned core with a 7. 62-centimeter BeO reflector indicated that even in the reflector 
the spectrum was still quite fast (fig. 12). Thermalized neutrons were essentially non- 
existent; even near the outer edge of the reflector the calculated thermal- to total -flux 
ratio (%/<Ptot) was 5xi0 " 4 - Mean neutron flux energy varied from 0. 41 MeV at the core 
center to 0. 07 MeV at the outer edge of the reflector, and thus indicated some softening 
of the spectrum. Spectral data are summarized in table VII. The median fission energy 

in the core was calculated to be 0. 35 MeV- Data on the spatial distribution of the flux, 

13 

summarized in table vm, indicate a leakage flux of 4x10 neutrons per square centi- 
meter per second. 
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Energy, eV 


figure 12. - Rux spectra in cylindrical, fast-spectrum reactor with 7. 62’Cent! meter BeO 
radial reflector* 4AEI spectra at core midpiane * ) 


Neutron gen eration time. - A neutron generation time, or prompt -neutron lifetime, 
of about 0. 3 microsecond was calculated for a representative power -flattened core using 
the a option in the DTF code. This technique is essentially the addition of an equivalent 
1/v absorber to the core followed by a k calculation. The neutron generation time 
for small Ak/k could then be calculated by 

Ak 


l 


P 



( 4 ) 


where 

U k -*v 
k ‘ kk v 

and 

k multiplication constant of core 

multiplication constant of core with 1/v absorber 
a. reciprocal reactor period, sec* 1 

No neutron lifetime data on a similar system were available for comparison to check 
the calculated value. However, small, fast, bare cores have neutron lifetimes of the 
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TABLE VE. - NEUTRON FLUX SPECTRA IN BeO REFLECTED FAST -SPECTRUM REACTOR a 


Energy range 

Core center 


Core edge 


Inside edge of reflector j 

Outside edge of reflector 


Fraction of 

Flics per j 

Fraction of 

Flux per 

Fraction of 

Flux per 

Fraction of 

Flux per 


neutrons 1 * 

unit 

neutrons 1 * 

unit 

neutrons 13 

unit 

neutrons* 3 

unit 




lethargy 



lethargy 


lethargy 


lethargy 

a. 7 to 15 MeV 

0,0421 

0.0301 

0,0341 

0.0244 

0, 0222 

0,0159 

0.0234 

0.0167 

0.82 to 3.7 MeV 

.282 

,188 

.234 

, 156 

. 176 

, 117 

,136 

,0907 

0.3 to 0.82 MeV 

.297 

.297 

.256 

.256 

.225 

.225 

.123 

.123 

0,11 to 0.3 MeV 

, 235 

.236 

.223 

.223 

.221 

.221 

.154 

.154 

0.025 to 0.11 MeV 

.125 

.0833 

,180 

.120 

.222 

.146 

.242 

,161 

5,5 to 25 keV 

,0166 

,0107 

,0608 

.0405 

.0972 

.0648 

.180 

, 120 

1.2 to 5. 5 keV 

9.2X10" 4 

6.1X10"* 

,0105 

.0070 

.0261 

. 0174 

, 0860 

; .0573 

0.28 to 1.2 keV 

7.7X10" 6 

5.1X10' 6 

.0011 

7. 3X10 " 4 

.0069 

.0046 

.0317 

.0211 

0.1 to 0.28 keV 

4.6x10"® 

4. 6x10"® 

1.4X10' 4 

1.4 


, 0017 

.0017 

.0100 

.0100 

0.03 to 0.1 keV 

Small 

Small 

2.0 


1,7 


7.9XKT 4 

e.exio" 4 

, 0062 

.0043 

11 to 30 eV 





Small 

Small ; 

3,9 

3.9 

.0023 

.0023 

3 to 11 eV 









2.6 

2,1 

,0017 

.0013 

0.88 to 3 eV 









3.7 

3.0 

.0015 

.0012 

0.41 to 0.88 eV 









2.5 

3, 3x10 ^ 

4. 3xi0“ 4 

5.7X10" 4 

0.1 to 0.41 eV 









9,5xi0~ 6 

6.7X10’ 6 

1.5 

1.1 

Thermal (<B. 1) 









5,9xi0" 5 


4.8 


Median energy, 


0.41 



0. 32 


0. 22 

0.075 

MeV 













■ 


^Calculated at midplane of 32,6-cm-dIam zoned core surrounded radially by 0.9 5- cm tungsten pressure vessel at 
7. 62 -cm BeO reflector. 

h r 15 MeV 

^Norraalifted to / qu(E)dE - 1. 


TABLE VHI. - SPATIAL FLUX DISTRIBUTION IN BeO 
REFLECTED FA ST -SPECTRUM REACTOR 


[Zoned core; 7.62-cm BeO radial reflector; 1.27 -cm W 
end plate; 0.95-cm W pressure vessel,] 


Fraction of neutrons absorbed 

0.67 

Fraction of neutrons leaking axially 

0,13 

Fraction of neutrons leaking radially 

0.20 

Flux at core center, neutrons/(cm^)(sec) 

8X10 13 

Flux at core edge (radial), neutrons/fcm^Hsec) 

4X10 13 
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TABLE EX. - TYPICAL REACTOR CHARACTERISTICS a 


Power, MW 

1 

Core life, hr 

25 000 

Coolant 

He- 60 Xe 

Coolant outlet temperature, K 

1650 

Coolant pressure, MN/m^ 

1,71 

Structural materials 

W and W-25Re 

Fuel material 

U 235 0 2 or U 235 N 

Radial reflectors, cm of BeO 

7,62 

Axial reflectors (end plates) , cm of W 

1.27 

Radial peak- to aver age -power ratio, 

1.05, 1.10, or 1, 35 

Axial peak- to average-power ratio, 

1.31 

Reactivity control range (unshielded), percent 

>10 

Maximum multiplication constant 

1.05 

Reactivity shutdown margin, percent 

>5 

Mean fission energy, MeV 

0.35 

Prompt -neutron lifetime, jisec 

0,3 

Flux core center, neutrons/(cm^){sec) 

8X10 13 

Leakage flux, neutrons/(cm^){sec) 

4xl0 13 


a Right- cylinder reactor core with length to diameter ratio of 1,0. 


order of 0. 02 microsecond (ref. 19), and small fast cores with thick Be reflectors have 
neutron lifetimes of about 5 microseconds (ref. 20). Thus, the calculated value of 0. 3 
microsecond for a fast core with an intermediate BeO-reflector thickness appears 
reasonable. 


Typical Reactor Properties 

The results of this study have indicated that a typical reactor design would have the 
properties summarized in table IX. 


SUMMARY OF RESULTS 

A model was conceived for a gas- cooled, compact nuclear reactor required to op- 
erate at 1 megawatt for 25 000 hours with a coolant outlet temperature of 1650 K. This 
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model has the following characteristics: 

(1) A cylindrical core with a length to diameter ratio of 1 surrounded by a 0. 95- 
centimeter -thick pressure vessel, axial end plates, and a radial reflector 

(2) A control system based on neutron leakage and regulated by axial movement erf 
a radial reflector 

(3) The materials or for fuel, W and W -25 Re for structural compon- 

ents, and BeO for reflector 

An analytical study of this reactor produced the following information: 

1. A nominal control range of 10 percent reactivity was required based on (a) cal- 
culated reactivity values of 1 percent for fuel burnup and 1. 8 percent for temperature 
effect, Os) an estimated 2. 2 percent reactivity for contingency, (c) an assumed zero 
Doppler coefficient, and (d) 5 percent reactivity for shutdown margin. Thus, the multi- 
plication constant for the reactor in the fully reflected condition is 1 . 05 . 

2. Reactivity worth of the radial reflector was shown to be a function of core diam- 
eter, reflector thickness, reflector material, and radial power distribution. For 
uniformly fueled, fully enriched cores with diameters less than or equal to 38 centi- 
meters a BeO reflector greater than or equal to 7. 6 centimeters does provide sufficient 
reactivity control range. The reflector worth for a 32. 6-centimeter-diameter core with 
a 7. 6 2 -centimeter -thick reflector increased by a factor of 1. 5 when the radial power dis- 
tribution was flattened by fuel enrichment variation to a peak- to average -power ratio of 
1. 05. Reflector worths were decreased to relative values of about 0. 8 when a 7. 62- 
centimeter BeO reflector was replaced by the nonmoderating materials Nb or Mo. 

3. A total BeO radial reflector thickness of 7. 62 centimeters was selected for this 
study. The reflector was divided into two sections, the inner section for shutdown con- 
trol and an outer section for operating control. 

4. For the selected end -plate thickness of 1. 27 centimeters, the calculated axial 
power distribution could be approximated by the function cos rrZ/l. 27 L which has a 
peak- to average -power ratio of 1.31. 

5. Calculations indicate that sufficient reactivity control range with a 7. 62- 
centimeter BeO radial reflector is available when a 4ir tungsten shield is placed adjacent 
to the reflector of a core zoned to a {P p /P) R of 1. 05. However, the effects on core 
design were not determined. 

6. Based on a typical core with a diameter of 32. 6 centimeters the following core 

characteristics were calculated: a median fission energy in the core of 0. 35 MeV, a 

13 

prompt -neutron lifetime of 0. 3 microsecond, and flux values of 8x10 neutrons per 

13 

square centimeter per second at the core center and 4x10 neutrons per square centi- 
meter per second at the core edge. 

7. The parametric data are combined into design curves relating core diameter to 
required fuel density. Three curves, one each for the power -flattened condition of 
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radial peak- to average -power ratios of 1. 05, 1.1, and 1. 35 (unzoned), were constructed 
for a cylindrical core with a 0. 95 -centimeter W-25Re pressure vessel, 1. 27 -centimeter 
W end plates, and a 7. 62-centimeter BeO radial reflector. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 26, 1968, 

120-27-06-05-22. 
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APPENDIX A 


CROSS SECTIONS 

Of the multigroup reactor codes available, the transport codes {such as the one- 
dimensional DTF and the two-dimensional TDSN) are particularly adaptable to small 
fast-spectrum cores. Also, the running time of DTF is sufficiently fast to permit 
parametric studies. However, before reactor calculations can be made it is necessary 
to obtain a set of neutron cross sections for the materials of interest. To accomplish 
this the GAM-n code was used to calculate P Q and Pj^ microscopic cross sections for 
14 neutron energy groups from 15 MeV to 0. 414 eV averaged over a fast flux spec- 
trum. The spectrum was generated by GAM -II in a medium composed of a nominal core 
composition, and the group energy split was similar to the Hansen -Roach cross-section 
set (ref. 21). Because of the importance of inelastic scattering in fast cores, the full- 
down scattering matrix in the GAM output was retained in the cross-section set. Also, a 
resonance calculation was performed to refine the spectrum in the eV to keV range. The 
energy range below 0. 414 eV was divided into two groups, and the cross sections were 
either generated by TEMPEST -II or hand calculated. No effort was made to refine these 
relatively rough values for the two low-energy regions because of their relative unimpor- 
tance to fast-reactor calculations. In addition, a set of Be cross sections was generated 
which were flux weighted over a moderating medium (BeO) spectrum. These cross 
sections were used in the reactor reflector regions. 

An initial check on the validity of the cross section was made by calculating a num- 
ber of critical assemblies (table X). These assemblies were selected because of their 
size and the similarity of materials to those expected to be used in this study. All cases 
were calculated with Pq transport -corrected cross sections, and the results were suffi- 
ciently close to k = 1. 0 to warrant confidence in their accuracy. The fact that all values 
(except the 710 critical) exceeded 1 is characteristic of the use of P Q transport- 
corrected cross sections. The low value of the 710 calculation probably is the result of 
a discrepancy in the transverse leakage. 

A more comprehensive check on the cross sections and the use of DTF was obtained 
by comparison of flux and power distribution .calculations with 710 critical experiment 
data. (Although specific data on this experiment are classified, it is similar to the 
model described in reference 20. ) Values dependent on a fast spectrum (central core 
region) were accurate, whereas those near the core edge (which were more dependent on 
low-energy neutrons) showed more deviation. At the core center the fission ratio of 
U 23 * 3 to U 235 was underpredicted by 2. 3 percent and the local- to minimum-power ratio 
was overpredicted by 10 percent, whereas at the core edge the local- to minimum -power 
ratio was under predicted by 29 percent. 
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TABLE X, * CRITICAL EXPERIMENT CALCULATIONS 


1 Core 

Calculated 
multiplication 
constant, , 

k 

Godiva 

1.012 

Be -reflected 
sphere 

1.042 

W -diluted bare 
sphere 

1.000 

W- diluted bare 
sphere 

1,014 

Ta -diluted bare 
sphere 

1.029 

W-reflected 

sphere 

1.009 

710 critical 
experiment 

.989 


Comments 


Small bare sphere of §3,9 percent enriched U 
(ref, 21) 

Small sphere of 93, 5 percent enriched U 
surrounded by a 5, 08 -cm Be shell 
(ref, 21) 

Bare sphere composed of 0,309 enriched U 
and 0.691 W (rel, 22) 

Bare sphere composed of 0,5 enriched U 
and 0, 5 W (ref, 23} 

Bare sphere composed of 0,495 enriched U 
and 0. 505 Ta (ref, 23) 

Small sphere of 93, 5 percent enriched U 
surrounded by a 10, 16 -cm W shell 
(ref. 21) 

Cylinder composed of 0,27 enriched U, 

0.26 W, 0, 1 Ta, and 0, 1 A1 with a 15-cm 
Be radial reflector (transverse leakage 
was estimated from a 2-D calculation) 

(ref, 20) 


However, since overall reactor properties in these small, fast cores are primarily 
dependent on fast spectrum properties, both the ealculational method and the cross 
sections appeared adequate for this study. An interesting aspect of the critical experi- 
ment data was the spectrum in the Be reflector region as defined by cadmium (Cd) ratios 
from foil experiments. Criticality calculations indicated that both the multiplication 
constant of the system and the spectrum in the Be region could be varied by separating 
the Be into two regions, one in which fast -spectrum -averaged cross-sections were used 
and one in which moderating spectrum cross sections were used. 

This division would seem logical since the fast spectrum from the core is slowed 
down (moderated) considerably by the Be, the degree of moderation being dependent on 
the depth of penetration. The use of Be fast cross sections over the entire Be region 
resulted in a larger multiplication constant and a lower Cd ratio. A calculation using Be 
fast cross sections for the first 4. 14 centimeters and Be moderated cross sections for 
the remainder of the Be region gave a closer approximation to the experimentally deter- 
mined Cd ratios. However, these were slightly high and it was decided to use Be fast 
cross sections over the first 5. 08 centimeters of the reflector region (BeO) and Be 
moderated cross sections for the remainder for the parametric calculations in this 
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design study, 
in table XL 


TABLE XI. - DEPENDENCE OF DTP CALCULATIONS ON 
INPUT OPTION 


Input option comparison 

Reactivity change, a 
percent reactivity 

s 4 -s a 

p 0 - P 1 4 

0. 15 
1* 4 

Convergence criterion, case with 1x10 
case with 1X10^ 

O 

Be cross sections, case with moderated 0 - 

-1*5 

case with fast 0 


Cross sections, case with resonance corrected - 

1*2 

case with nonresonance d 



a Based on criticality calculations of 710 critical experiment. 

^Transport -corrected Pq cross sections. 

c Type of spectrum weighting in GAM -II calculation. 

^In a nonresonance calculation GAM-n uses a resonance integral for 
calculating cross sections. 


The effect of Be cross sections on the multiplication constant is indicated 


27 


I 





APPENDIX B 


CALCULATION TECHNIQUES 

A few of the options in the use of DTF were investigated prior to undertaking param- 
etric criticality calculations (table XI). Based on these results all reactor calculations 
were made using the transport -corrected PqS 4 code option with a convergence criteria 
of lxio'4 and resonance -corrected cross sections. Although this procedure indicates 
that multiplication constants could be overpredicted by i. 55 percent Ak, the low value 
calculated for the 710 critical experiment (-1.1 percent Ak in table X) represents a 
counteracting correction. Therefore, multiplication constants were reported in this 
study as calculated. 

The calculation model used homogenized regions in cylindrical geometry and, for 
one -dimensional problems, a buckling to simulate transverse leakage. The buckling is 
calculated internally in DTF by using an extrapolated dimension; for ©cample, for axial 
leakage from a cylinder this equals the actual height plus reflector savings plus twice an 
extrapolation length of 0. 7l/£ tr . 
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